QUANTUM PROBABILITY MEASURES AND 
TOMOGRAPHIC PROBABILITY DENSITIES 



G.G. Amosov^ and V.I. Man'ko ' 

Abstract. Introduced recently approach based on tomographic 
probability distribution of quantum states is shown to be closely 
related with the known notion of the quantum probability mea- 
sures discussed in quantum information theory and positive oper- 
ator valued measures approach. Partial derivative of the distribu- 
tion function of quantum probability measure associated with the 
homodyne quadrature (symplectic quantum measure) is shown to 
be equal the tomogram of the quantum state. Analogous relation 
of the spin tomogram to quantum probability measure associated 
with spin state is obtained. Star-product of symplectic quantum 
measures is studied. Evolution equation for symplectic quantum 
measures is derived. 

1. Introduction 

Recently [H 12 E] the tomographic probabihty distributions were 
found to be related to Wigner functions The tomoraphic probabil- 
ity distribution of E] was used in optical tomography scheme E] 
to reconstruct the Wigner function of photon states by measuring the 
homodyne quadrature distributions and applying the Radon transform 
[7] to find the Wigner function. In |3j an extension of optical tomogra- 
phy to sjTiiplectic tomography [S] was suggested. In framework of the 
symplectic tomography scheme the Wigner function and density opera- 
tor can be reconstructed using Fourier - like integral (instead of Radon 
integral) of the symplectic tomogram. In [0] the tomographic approach 
was shown to be connected with well-known star-product quantization 
procedure f[TU]). 

In quantum information theory and positive operator valued mea- 
sures approach (in the context of quantum measurements) [llj the 
notion of quantum probability measure for a generic variable is a basic 
concept. Namely, to each pair (p, M) consisting of the quantum state 
(the density operator) p and the positive operator valued measure M, 
which can be identifyied with the quantum observable, it is associated 
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axiomatically (in the way of ^2]) the probabihty measure /i^^ which 
determines the readings of a classical measuring instrument. Taking 
into account the enssemble {/^^} for the fixed measure M and all states 

p, it is possible to reconstruct M. In quantum information theory 
the states p play a role of quantum probability measures used to en- 
code the information. After the states p are transmitted through the 
quantum channel, one uses M as desicion rules allowing to decode the 
data from the channel output. It is worthy to point out that one can 
connect a state with a probability distribution density for arbitrary ob- 
servable a using the characteristic function techniques (see jE]). This 
construction is related to the construction given by von Neumann in 
jl5j . which is used in quantum information theory pTl I13j . 

The aim of our work is to establish a new connection of tomographic 
probability with such well studied mathematical object as quantum 
probability theory in addition to known connection with the star- 
product. Namely, we show that the symplectic tomogram can be ob- 
tained as a partial derivative of the distribution function of quantum 
probability measure. Till now the explicit relations between the to- 
mographic probability distribution (tomogram) of quantum states and 
quantum probability measure concept (for the specific observable) used 
in quantum information theory were not known. 

The paper is organized as follows: In Section 2 we review the quan- 
tum probability theory. In Section 3 and 4 symplectic tomography and 
spin tomography respectively will be discussed. In Section 5 connec- 
tion of tomogram with quantum probability measure will be obtained. 
In Section 6 new concept of star-product for symplectic quantum mea- 
sures will be introduced. In Section 7 evolution equation for symplectic 
quantum measures will be obtained. Srction 8 is devoted to conclusions 
and perspectives. 

2. Basic concepts of quantum probability 

Let if be a separable Hilbert space. A positive linear operator p in 
H is said to be a state (density operator) if Trp = 1. We denote C{H), 
C+{H) C C{H), <j{H) and / the sets of hermitian opeartors (quantum 
observables) , positive operators, states and the identity operator in H 
correspondingly. A map M transmitting each Borel subset f2 C M to 
a positive operator M(fi) G C+{H) is said to be a positive operator 
valued measure (POVM) if ^ M{Qi) = I for any finite or countable 

i 

fragmentation M = Ufij in the sense of the strong operator convergence 
of the series. The POVM M is said to be orthogonal if M(fi)2 = M{n). 
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By means of the spectral theorem, given a G C{H) there exists an 
orthogonal POVM M such that 

(1) a= [ XdM{{-oo,X]) = [ X5{d- X)dX, 



where S{d — X) is a density of the measure M. It means that the 
measure M is expressed in terms of the Heaviside function 

M((-oo, X]) = e{d-X) = ^ [ le^"'^^-''^dk. 

2m J k 

The presentation of the spectral theorem by means of the operator - 
valued delta function can be understood as one can check that matrix 
elements of last two parts of (Q) in eigenvector basis of the operator d 
coincide if we use standard Dirac delta function properties. In this way, 
it is possible to identify observables d G C{H) with orthogonal POVM. 
We shall say that d G JC{H) is associated with the POVM M appearing 
in its spectral representation. Taking p G (t{H) and the POVM M one 
can put 

(2) Mf{n) = TripMm 

for any Borel set C M. Then Ai^^ is a classical probability measure 
on R. Calculating the matrix p{x, x') of the density operator p in the 
basis in which d is diagonal, we can represent Q in the form 

(3) Mfi^l) = I p{x,x)dx. 



If the POVM M is associated with some observable d G jC{H), we shall 
write Ai'i as well as Ai'-^ . In quantum probability theory ^T], it is 

postulated that Ai^ = defined by the formula (j21) gives us a distri- 
bution of the readings of a classical instrument measuring the quantum 
observable d associated with POVM M in the state p. Notice that for 
the orthogonal POVM the approach introduced above appeared firstly 
in jT3] . The same probability distribution can be obtained by means of 
the characteristic functions techniques J3]. Alternatively, in quantum 
information theory [12] M determines certain decision rule allowing to 
decode the information containing in the state p. 

For example, consider the position operator x G C{H) acting in 
the Hilbert space H = L^(M) by the fromula {xf){x) = xf{x), / G H. 
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Denote xn the characteristic (indicator) function of the Borel set f2 C M 
such that 

Involve the operator xn acting by the formula {xnf){x) = 
Xn{x)f{x), f E H. Then the map M defined for a Borel set Q by 
the condition M{Q) = x^ is an orthogonal POVM on R. One can see 
that M gives us the spectral decomposition of the operator x such that 



(4) X = J Xdx(^oo,x]- 

R 

Denote \X >< X\ the density of the orthogonal POVM x, then 
transforms to the following, 

x = [ X\X >< X\dX, 



\X >< X\dX = i. 



In this way, the probability measure determined by the state p can be 
written as 



(5) Mj{n) = Tr{pxn) = j p{X, X)dX, 

where p{X, X') is a matrix of the density operator p in the basis con- 
sisting of generalized eienvectors \X >< X\. For a pure vacuum state 
Po = iV'o >< ''/'ol in the Schrodinger representation 

1 

the formula (0) yields 

(6) M%{n) = exp{~x')dx, 

n 

which is a Gaussian measure with the zero mean and the variance equal 
to 1/2. 

f - \ 

Let Jz ~ { Q 1 1 be an operator of spin projection on the 2;-axis 
for the particle with the total spin J = h- Picking up the Euler angles 
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ip, 9 one can define the rotation matrix by the formula 



(7) R{<p.i^.e)= .''"'f 



cosi^e 2 istn-e 2 
tstn^e 2 cos|e 2 



Notice that ((7j) determines irreducible representation of the group 
SU{2). It is straihtforward to check that a spectral decomposition 
of the operator a = R{(t), ip, 6)JzR{(t), ip, 0)~^ is given by the formula 

\sinee'^ -\cose J ~ 2\ ^sinOe^^ sin^^ 
1 f sin^l yinOe-'^ 



2 V -^sinOe'^ cos^^ 



Using (jS)) we obtain for the probability measure which is discrete 
in the state p = \ip >< "ipl with li^ >= ^ g ^ the distribution function 



0, x<-i 



(9) Fix) = M%{-oo,x]) = { cos% -\<x<\ 

The function Q gives us the Bernoulli distribution concentrated in 
two points Xi = — I and X2 = \ with the probabilities p = cos^^ and 
l—p = sin^^ correspondingly. 



3. Tomographic representations for continuous variables 

Tomographic probability (tomogram) determining a quantum state 
is introduced by relation 

(10) ^{X, /i, z/) =< 6{X — fix — up) 

where x and p are position and momentum operators and the density 
operator p defines the averaging for arbitrary observable a by 

(11) < a >p= Tr{pa). 

The Dirac (5-function in (fTUI) is defined by its Fourier decomposition as 

' Aka , 



(12) 6{a) = — / e'^^dk. 

2vr J 

According to [SI El E] the tomogram ()10|) is related to the density matrix 
in position representation piy, y') by 

Iff i(y'^-y''^)fi iX(y-y') 

(13) w{X,p,u) = ^— p{y,y)e 2. dydy 
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and with the Wigner function defined in |3] as 

/?/ 11 

via the relation 

(15) w{X,f^,u) = J Jw{q,p)5{X-f,q-up)^. 
The tomogram determines the Wigner function as 

(16) Wiq,p) = ^j j j w{X,fx,u)e'^^'f"'-''P^dXdfxdu. 

It means that the tomogram can be used to describe the quantum states 
completely. In terms of a wave function ip{x) the tomogram reads for 
the pure state p = {ip >< [H 



(17) w{X,fi,u) = -^1 / ij{y)e^-'^dy\' 



Put 



Pn = \lpn >< 1pn\, U = 0, 1,2, 



(18) < A-|V'„ >= ^-^i/„(A-)e.p(-^) 

are wave functions of the excited state of an oscillator. Here Hn, n = 
0, 1, 2, . . . are the Hermite polynomials. Then (fT7j) gives us the tomo- 
gram for the state /3„ as 

111 X X"^ 

Let us involve the wave function of a coherent state iV'a >< ^'al; a G C, 

< X\iPa >= J^'^'^Pi-^ + - Y - 

Then the tomogram of the state {ipa >< i^al given by 



1 f - V2Reafi - V2I 



(20) ./;«(X) = ^ ^ y.. - v....^p- v..„.^.)^ 



The most important property of the tomogram is that it is a standard 
density of the probability distribution function on M, i.e. 

(21) w{X,fx,iy)>0 
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and 

(22) / w{X,fi,u)dX = 1. 



The physical meaning of the real variable X is that this variable is 
equal to position of a particle measured in rotated and scaled reference 
frame of the particle phase space, in which we get 

(23) X = jiq + up, jj, = e^coscj), v = e'^sincj), 

where and A are an angle of the rotation and a real sqeezing parameter 
corresp ondingly. 

4. Tomography of spin 

Here we review the tomogram of discrete variables (spin). The to- 
mogram of a spin state is defined by the relation 

(24) w'-^\m,e,(f)) =< 6{m ~ R^^\(f),'^,e)J,R'-^\(f),^,e)-^) >p, 

where S is the Kronecker operator delta function which for arbitrary 
hermitian operator u with integer eigenvalues is given by the Fourier 
integral 

27r 





a nonnegative half-integer j is the total spin, m = —j, + is 
the spin projection on ^-axis, is a operator of spin projection on the 
z-axis and R^^\(j), if), 9) is the matrix of irreducible representation of the 
group SU{2) in standard basis, elements of the group are parametrized 
by the Euler angles (f),ilj,6. The state of spin is determined by the 
hermitian nonnegative 2j + 1 x 2j + 1 - density matrix with unit trace. 
One can show ^Tj that the density operator can be found from 
fl24|) . In fact, the formula (j^^ can be considered as a linear system of 
equations determining the density matrix if the tomogram is a known 
function. 

5. Quantum probability and tomograms 

By the definition of the orthogonal POVM M given by the spectral 
decomposition of the operator d, we get -^M{{—oo,X]) = S{X — a). 

Suppose that M is determined by the spectral decomposition of fix+up. 
Comparing dTUl-fTTl) and we get 

(25) wiX,^,,l.) = ±Mf ^"^ii-ocX]). 
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Using the relation ()25|1 one can calculate the distribution function of 
the measure M'f''^'^^ for the states pn (fTHj) in which the tomograms Wn 
are given by the formula (fTIH) such that 



X ^ 
(26) , ^ I Hli , ^ )exp( J^)dy. 



Take the function 



2 a\ 



1/4 2 



and decompose it into the series over the excited wave functions of an 
oscillator ipn{x) (fT5|) . then 



n=0 

It follows that 

X X 



7~Tn ZTn V n\ V ml J J 



n=0 m=0 

X 



-OD 



^ y exp(-?/ + V2{a + (3)y )dy 



(27) V(l + er/(X-^)). 

Using we obtain for the probability measure (j^Ull 

A^^f+^^((-oo,X]) = 



In particular, for the vacuum state po we have 

A^J+-((-oo, X]) = i(l + er/( ^^ )), 

which is in accordance with (jH]) if one put /i = 1, z/ = 0. 
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In this way, using the tomogram it is possible to obtain the proba- 
bihty measure Ai^^^^^ without a calculation of the spectral decompo- 
sition for the operator /xx + vp. We shall call Ai^^^^'^^ by a sympectic 
quantum probability measure. 

Analogously, suppose that the spectral decomposition of the operator 
a = e)J-,R^i\<i), ^/j, 6)-^ has the following form 

j 

k=-j 

where M^, —j<k< j, form the orthogonal POVM. It follows from 
(1^ that 

j 

k=-j 

J2 ^kmM%{k - 1, k]) = M%{m - 1, m]), 

k=-j 

where 6km is the Kronecker symbol, m = —j, + 



6. Star-product for symplectic quantum measures 

To formulate quantum mechanics using a map a fa{S) between 
the set of hermitian operators and the set of functions on appropriate 
space, one needs to involve a new multiplication rule for the functions 
which should be associative. This multiplication is said to be a star- 
product (see !T0^ ) and, in particular, it can be involved for the functions 
defined by the map (x = (X, yU, z^)) 

(28) fa{x) = Tr{a5{X -fix- up)) 
in the way of 9J such that 

(29) fa * fl{x) = j K{xi,X2,x)fa{Xl)fl{x2)dXidX2, 

where the kernel K{xi, X2, x) is written as 

K(xi,X2, x) = Tr{V(xi)'D{x2)S{X — /ix — up)) 



and dx = dXdiidv. Here the operator V{x) is defined by the formula 
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The operator d can be reconstructed from the function fa{x) by the 
formula 

(30) d = J J J U{x)V{x)dXdfidu. 

In particular, if d = p G cr{H) is a state, then ()28|1 yields a number of 
the probability distribution densities associated with p. 
Now consider the map 

(31) p^Mf,ix) = Mf'-'^i-oo,X]) 

from the set of states (y{H) to the set of probability measures on the real 
line. The spectral theorem determines a representation of arbitrary ob- 
servable d as an integral over pure states. Hence it is possible to extend 
the map (jHTj) from the states p to all observables d by means of a linear- 
ity. Taking the spectral decomposition px + up = J XdM{{—oo, X]) 
we get 

(32) Ma{x) = Tr(dM((-oo,X]))). 

Here for an arbitrary observable d we obtain the family Aia{x) deter- 
mining a non-positive Borel measure on the real line for any fixed p and 
u. We shall call them by symplectic quantum measures. Comparing 
(|Sni),(I23) and dng) we obtain 

(33) = ^'^(^^ 

and the operator d can be reconstructed from the family of symplectic 
quantum measures Aia such that 



(34) d = / / / V{x)dMa{X)dpdu. 




Here dAia{X) denotes the measure on the real line determined by 
Aid{x) for the fixed p and z/. In the formula (jS^ and below we claim 
that, at first, the integration is done over dA4d{X) for the fixed p, v and 
then we integrate over dp and dv. The formula (j34p allows to define 
an associative multiplication for the symplectic quantum measures M.a 
which can be derived from (P^ . 

(35) M-a-^Mi{x) = j k{xi,X2,x)dMdiXi)dpiduidMi{X2)dp2du2, 

where the kernel K{xi,X2,x) defined by the formula 

K{x,,X2, x) = Tr(P(fi)P(f2)M((-oo, X])), 
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xi = (Xi, /ii, z/i), X2 = {X2, fi2,^2)- By a construction the multiphca- 
tion (j35p is in accordance with the obvious muhiphcation on the set of 
observables, such that 

(36) M,f,{x)=Ma*M-,ix). 

It follows from ()36p that the symplectic probability measures Aip asso- 
ciated with the pure states p = \ip >< ipl are idempotents with respect 
to the multiplication *, i.e. 



7. Evolution for symplectic quantum measures 

It was shown in 19^ that if the quantum state satisfies the time- 
evolution equation 

(37) dtp=-t[H,p], 



^2 



then for systems with Hamiltonian of the form H = y + V"(g) the 
evolution of the quantum tomogram (fTO|l by virtue of (|37jl is given as 

dw dw 
dt ^ dv 



(38) -,(^(-(_)----_)-n-(^)-5^ + ^^)K- = o. 

It is straightforward to check that the equation (jHHj) determines the 
evolution of the functions (j^Hj) and, therefore, by means of the connec- 
tion of the distribution functions Aia{X) of symplectic quantum 
measures (j32|) . Integrating (j38|) over X we obtain the following form 
of evolution, 

dMajX) dMgjX) 

dt ^ du 

(39) 

-'-<-<^'-'| - f ^' - ^'H^'-'l + f ^)'^« 

Note that we have derived the linear evolution equation valid only for 
symplectic quantum measures. The general equation for continuous 
nondemolition measurement introduced in [201 E] is nonlinear. 
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8. Conclusion 

To conclude we summarise the main results of the paper. We estab- 
lished relation (see Equ. ((23)) of symplectic tomograms to distribution 
functions of symplectic quantum probability measures, which gives a 
possibility to use results of the quantum probability theory to study 
properties of quantum tomograms. On the base of the established rela- 
tion the evolution equation for the symplectic (non-positive in general) 
quantum measures is obtained (see Equ. (jHUjl ). The star product of 
symplectic quantum measures is introduced (see Equ. (jSSl))- The ap- 
proach is illustrated by example of several different states of harmonic 
oscillator. In addition, the quantum probability measures and its con- 
nection with tomogram are constructed for such discrete observable as 
spin. 
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